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I develop here the microscopic quantum theory for description of creation of phonons and rotons 
in superfluid helium by a solid heater. Starting with correct transfer Hamiltonian describing a 
coupling between the solid and liquid 4 He the probabilities of transformation of a single phonon in 
the solid into i) single helium phonon, ii) two helium phonons, and iii) single helium roton are found 
out. All the obtained expressions account for different polarizations of phonons in the solid. The 
heat transfer associated with single phonon and single roton channels are calculated. Particularly, 
the obtained expression for heat flux via the single phonon channel calculated in the framework 
of present microscopic theory exactly coincides with the well known Khalatnikov formula obtained 
initially in the framework of acoustic-mismatch theory. The impossibility of direct creation of Pj - ' 
rotons becomes clear in the used framework due to accurate account of the boundary conditions at 
the solid - liquid helium interface, which is in agreement with recent experimental results. 
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I. INTRODUCTION 

The elementary excitations in superfluid 4 He have been the object of extensive study for many years. The phonon- 
roton type excitation spectrum proposed phenomenologycally by Landaun allows to account correctly for superfluid 
properties and thermodynamic behaviour of liquid helium below the A-transition temperature. The existence of such a 
spectrum was demonstrated by the numerous direct observations of single excitations from inelastic neutron scattering. 
In spite of many various attempts to built a completely microscopic theory to derive directly the thermal excitations 
spectrum, the general picture drawn first by Landauu and Feynman and Cohcno had not changed essentially during 
more than four decades. 

Since we realized the general picture of elementary thermal excitations in bulk liquid helium, the natural problem 
appeared: how to describe a heat flow across a boundary between a solid and liquid helium. All earlier investigations 
of this process were concentrated on studying of a discontinuity in temperature appeared between the two materialsB, 
which is commonly known now as Kapitza resistance. The first consistent accurate theory of this effect was developed 
by Khalatnikovu. He built his theory starting with quantizing of elastic waves of different polarizations in the solid. 
But in general, the main results of Khalatnikov's theory are completely consistent with classical acoustic-mismatch 
theory due to LittleO. The in-depth reviews of theoretical and experimental investigations of Kapitza resistance 
problem are given by Challiscl and recently by Zinov'cvaQ. In their pure form all these theories show that transfer 
of phonons from the solid to the liquid helium is impeded first due to a large difference in densities pl/ps <C 1 (pl 
and ps is the densities of liquid helium and solid, respectively), and second due to the smallness of critical angle 
sm8 c = c/s t <C 1, where c is the first sound velocity in helium and St is the velocity of transverse sound in solid. The 
total internal reflection occurs when the angle of incidence exceeds this angle 6 C . Because of this the helium phonon 
can propagate only into a narrow cone along the direction perpendicular to the interface. 

According to the present understanding, solid can inject the phonons in two groupsa. The first is phonons in the 
narrow cone perpendicular to the solid heater surface, while the second group is phonons emitted into the entire half 
space, these are the so-calledr-background phonons. The existence of the first group is understood satisfactory from 
acoustic-mismatch theoriescTQ. This mechanism was directly confirmed by experiments^ (see also Ref. 7). Howayer 
the existence of background phonons, which is evident from experiments on phonon-emission from cleaved crystalsoEEl, 
was not satisfactory explained theoretically until now in spite of many attempts! 1 ]. In fact, all these theories were 
devoted to various modifications of the acoustic-mismatch theory to account the imperfection of the interface and 
some additional scattering mechanisms. As a result, only small corrections to the Khalatnikov theory were obtained. 

All theoretical works cited above have been based on the classical acoustic theory and there was thus need for a 
microscopic quantum mechanical description since 4 He is a quantum liquid rather than a classical elastic fluid. There 
were several attempts to build a quantum mechanical theory of solid-liquid helium heat transmission based on various 
simplified, modelalj, but the first really essential contribution to the subject has been done by Sheard, Bowley, and 
Toombsliij (hereinafter is referred as SBT). In fact I present here the modification of this theory. 
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The SBT theory is based on the description of coupling between the solid and the liquid 4 He via a transfer 
Hamiltonian in analogy with the theory of electron tunneling through a metal-insulator-metal junction. The main 
insight of Ref. 13 consist in the expression of tunnel Hamiltonian through the displacement of the solid at the interface 
and the stress tensor of liquid helium. 

The main shortcoming of SBT theory is the clear internal asymmetry of their transfer Hamiltonian with respect to 
solid and He II. Namely, as it follows from their theory, the solid acts on the helium, but there is no return action. 
As a result the transfer coefficient for the, single phonon channel (single solid phonon - single helium phonon) differs 
from the result obtained by KhalatnikovQ, which is correct for this process without any doubt. The authors of Ref. 
13 have also roughly estimated the transfer coefficients for single solid phonon - two helium phonons (two phonons 
process) and single solid phonon - single roton processes. As it will be shown below the more accurate computation 
of these coefficients gives rather different results. 

Calpilation of probability for creation of rotons by a solid heater is called for by recent experiments of Exeter 
groupli-3, which obtained quite a number of very interesting results. Particularly, they concluded that creation of 
the R( + ) rotons (rotons with wave vector greater than corresponding to roton minimum fco ~ 1.92 A) docs not 
essentially depend on the particular geometry of the heater, while the R^ - ) rotons (rotons with wave vector smaller 
than corresponding to roton minimum fco ~ 1.92A) cannot be created by a simple flat heater at all. This observation 
still requires the proper explanation. 

I this paper I propose a microscopic quantum theory for creation of thermal excitations in superfluid He by a 
heated solid. The starting point of the developed theory is the construction of more correct transfer Hamiltonian 
to exploit further the main idea of SBT theory!! To test the developed theory I first consider the heat flow via the 
single phonon channel with account of all phonon polarizations in a solid. The obtained expression for a heat flux 
through the interface associated with this process coincides exactly with the well known Khalatnikov resultcl. This 
direct support of the formalism allowed me to consider also the single roton emission process as well as to write down 
the probabilities for emission of a phonon pair by a single phonon in the solid. Particularly, I argue that emission 
of R(") rotons is forbidden, by the boundary conditions at the solid - He II interface, which is in agreement with 
experimental observationaiJ. 

The paper is constructed as follows. The transfer Hamiltonian corresponding to solid - liquid helium interface is 
constructed in Section II. In Section III I write down the tunnel Hamiltonians for the particular single excitation and 
double excitations processes. The well-known Khalatnikov formula for a transfer of single solid phonon into single 
phonon in helium is obtained from the general formalism in Section IV. Finally in Section V I study the creation of 
rotons from a heated solid via the single excitation channel. The outlines and conclusions are given in Section VI. 



II. SOLID-LIQUID COUPLING HAMILTONIAN 



The whole system under consideration consists of the solid and the superfluid helium with a common boundary 
between them. To study the excitations transmission processes through the interface we must pick up the part of 
total Hamiltonian, which corresponds to the interface itself. Just this interface Hamiltonian represents the coupling 
between two media. 

Let ^(r) and 4>{r) be the field operators (creation and annihilation operators in Scrodinger representation) in 
helium. Here r is the coordinate of Lagrange point in superfluid helium. Suppose that thermal perturbation lead to 
the displacement of some atom in solid with equilibrium coordinate R„ by vector u„. In the frame of pair interaction 
approximation we can write the deviation of a total Hamiltonian from the equilibrium one in the following form: 

Hint = f ^(r) #int (|R„ + u„ - r - 5r\) ^(r)d 3 r, (1) 

where J7l is the volume of the helium, ?7 int is the potential of interaction between a physical point of helium with 
a volume element in solid. This potential behaves like a single-particle operator when act on the creation and 
annihilation operators in helium. The quantity Sr in eq. (|l|) stands for the physical point displacement in helium due 
to the perturbation u„ in solid. It is connected with the density fluctuations in helium by the relation 

On 

5n(r) = n(r) - n a = -^-Sr. (2) 
or 

Let us assume the interface to be a plane perpendicular to zth direction as is shown in Fig. |l|. Supposing the 
displacements u and Sr to be small, we can rewrite expression (lit) up to the first order in these perturbations as 
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Wint = J d 3 rn(r)U (r) + J d 3 ru z (x,y,0)f( SL \r)~Jd 3 r6r z (x,y,0)f ( z LS \r), (3) 



Or Or, 



where 



f( SL )(r)=V-t (r) ^ VRn r>. nt(Rn „ r) ^ (r)j (4) 

71 

f< LS )(r) - v-t (r)Vr ^r>. nt(Rn „ r)v ; (r) (5) 

n 

are the forces per unit volume of liquid helium and solid, respectively, and 

U (r) =^Ui nt (R„-r). 

n 

In the expression ([^) we take the quantities u z (x,y,0) and Sr z (x,y, 0) at the plane z — Z = because of assumed 
short-range interaction. The first term in eq. (||) corresponds to a static interaction and does not depend on the 
local spatial perturbations. Thus, it naturally should be included to the equilibrium Hamiltonian. To transform 
remaining two terms we have to account correctly the boundary conditions at the interface. First, require the normal 
components of displacements in solid and liquid helium at the surface to be equal 

u g (x,y,Q) = 6r z (x,y,0). (6) 

Second, the total work done must be zero. This requirement can be expressed mathematically as 

J dzfi 3 ^ (x,y,z)+[ dZfW (X, Y,Z) = 0. (7) 

Using the relations (||) and ([?]) we can express the interaction Hamiltonian (||) in the following simple form 

7i int = 2 / d\u z {x, yi Q)fi SL \v). (8) 



Note, the interaction Hamiltonian (|^) obtained here differs from SBT Hamiltoniariij by the factor 2. This difference 
follows from incorrect account of mutual influence of solid and liquid helium by the authors of Ref. 13. They have 
accounted only the perturbation in the solid, which leads to an incident excitation. But as it can be seen from the 
expression ([!]) above, the account of corresponding perturbation in helium leads to a return response, which doubles 
the interaction Hamiltonian in view of the boundary conditions (0) and (Q) . 
Using the relation 

V R „J7 int (R„ - r) = -V r #tat(r - R„) 
and the expression (g) for f < si ) we can get the following relation 

/(^) (w )«-r>(r)^°, (9) 
az 

where n(r) is the number density of helium atoms. Now we come to the first essential approximation. Let us choose 
Uq(z) in the form of repulsive potential step at the solid-liquid interface 



DbW - { l°{ 



Vo (z<0) nn s 
z>0) ' [W) 



wheije-|Vb = const > 0. Below I will use the method of elimination of the force (g) from Hamiltonian (g) proposed by 
SBTtS Due to their insight the height Vq of the step ( ^0| ) is not crucial. Note, the combination —n(r)(dUo/dz) is the 
external force density exerted on the liquid helium by the solid. The hydrodynamic equation of motion looks 
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where j is the momentum density of liquid helium, irj-i is the fluid stress tensor. On the other hand the quantum 
mechanical equation for evolution of operator j (Heisenberg equation) gives 



dU <9tt 2/ 

n(r) ^7 + fe? (12) 



where Hl is the Hamiltonian of liquid helium. Following Ref. 13 let us integrate equation (|12|) in the interval 
z G {— e, +s}, where e is of the order of distance over which an average helium density decays from its bulk value to 
zero inside the potential step. In fact, I will put finally e — > 0, which is consistent with our choice of the potential 
(|Io|). Since the momentum density j remains finite the integral over the left hand side of eq. (|l2|) behaves like o(z) 
and therefore can be neglected. The derivatives of all components of vector ir z i except d-n zz /dz are finite as well. So, 
we get 

+e 

- I n(r)—^dz = tt zz -n zz = n zz (x, y, 0). (13) 

J aZ z—e z ——e 

— e 

Finally, using the formulae (||), (^|), and ( |l3|) we can write the interaction Hamiltonian in the form 

Hint = 2 / u z (x,y,0)ir zz (x,y,0)dxdy, (14) 



where £ is the area of solid-liquid interface. Formula ( |l4| ) represents the interaction Hamiltonian, expressed through 
the solid displacement at the surface and diagonal component of liquid stress tensor. Now it should be represented 
in terms of collective variables, which allow to carry out second quantization in the simplest way. It is naturally to 
present ir zz as a sum of kinetic and potential parts as follows 

7r„=#)+7r^. (15) 
The explicit forms oLjf z k z n ^ and ir^ ^ are derived in Ref. 13 using the semi-microscopic theory of superfluids developed 



by Sunakawa et. al£3. Thus, in our notations the kinetic part of (15) takes the form 

rkinl 2 h 2 f dn _ 1 dn d 2 n\ ^ h 2 k 2 

zz = mnv 2 z + — — n l —- — )- y ~w~7^~- (16) 



2m V dz dz dz 2 J ^ 2mO.L 

k 

where m is the mass of helium atom, n is the number density operator, which stands for the first collective variable, 
v z is zth component of the velocity operator, which stands for the second collective variable, k is the wave vector of 
excitation in liquid helium. 

The potential part of (^) can be obtained from the equation of motion, which leads to the relation 

E = { / d 3 r'^ (r')V(r - r')^(r') J V-(r), (17) 

k ^ ^ 

where V(r — r') is the inter-atomic potential in helium. In fact the equation ( jlTj) implies the pair interaction, which 
strictly speaking is inappropriate for the liquid helium. We however will eliminate it from all final relations by the 

use of semi-phenomenological notations associated with superfluid helium. The detailed derivation of ir^ ^ from eq. 
( |l7| ) can be found in Ref. 13. 

In terms of Fourier transforms defined by v = ^ k Vkexp(ikr) and n = ^ k exp(— ikr) the kinetic and potential 
parts of zth diagonal element of the stress tensor up to the second order in collective variables look 



k,k' 



h 2 % 2 h 2 
— Y^ n k e- ikr -V^-^, (18) 
4m ^ z ^ 2mft L ' v ; 

k k 



r (P0t) 



k. P 

p^O 



Vz 



(19) 
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where V(k) is the Fourier transform of interatomic interaction potential in helium, which will be expressed below 
in terms of excitation spectrum of superfluid helium. Above we have used only two essential approximations: the 
step shape of potential ( fTo| ) and the pair interaction between helium atoms which results in the expression ( |l9| ) . The 
first approximation can be naturally justified in view of smallness of boundary inhomogeneity domain, which is of 
the order of interatomic distancestS. To avoid the explicit use of the second approximation I will act as follows. The 
most general approach to the problem of low lying energy levels in superfluids lead to the following form of excitation 
spectrumlill 

h 2 k 2 

where Tilo^ is the energy of excitation with momentum fik, is the structure factor of superfluid helium. Within the 
framework of microscopic theorytLa using the approximation of pair interaction this structure factor takes the form 



h 2 k 2 J v ' 



This particularly means that if we know the correct excitation spectrum of superfluid helium, say from experiment, 
then the quantity V(k) can be unambiguously expressed through the parameters of this spectrum. Namely 

w-rffi -***)• (22) 

where the analytic form of Wk depends on particular range of wave vectors. In this manner the problem becomes 
self-consistent. 

The Fourier transforms of collective variables can now be written in the second quantized forma 



/ A \ ^ Tik f p \ ^ 

" k = 1 J (" k _ aLk ) = 1 m \?£tn L ) (" k " a - k ) ■ (23) 



v k = -i-———— (ak + at k ) = -if — ] rfak + a^k). (24) 

where a k and &k are the creation and annihilation operators of elementary excitation in liquid helium with momentum 
Kk and energy Wk- 

Now we must write down the zth component of lattice displacement in the solid in terms of creation 6^ and 
annihilation bq a operators for solid phonons. Here q Q is the wave vector of solid phonon with ath polarization. 
Taking into account the longitudinal and transverse polarizations of bulk phonons and surface phonons we can write 

«*<*> - (^) " 2 fe « (s - K) + E ^ («. - *,) 

E^ 1 ^-^)^ 11 }' (25) 



+ 



where ps and f2s are the mass density and volume of the solid, respectively, u qa is the frequency of ath phonon 
mode, q s is the two-dimensional wave vector of surface phonon, rii = R|| is the parallel component of vector r (or R) 
at z = Z = 0. The angular amplitudes B a {0) can be determined aso 



. , _ 2s 2 cosg cos29 t 

l[) s?sin2^sin2^o + s2 C os 2 2^' 



m = 2s 2 cos 0q sin 29[ 

tK ' s 2 sin 28 l sin 26 + s 2 cos 2 26 ' K ' 
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(28) 



_ Qa + 4 [ + «t K * 9s + K t ^gs 



Wo 



Vc- 2 - s« 2 , (a = l,t). (30) 



Here c is the velocity of first sound in superfluid helium, si and s* are the velocities of longitudinal and transverse 
sounds in the solid, 6$ and 6 a arc the incidence and reflection angles. 

In fact, the formulae (|l4|), (p^)-(|2l]), and ( |25| ) give us the sought transfer Hamiltonian in second quantized form, 
which describes the transmission of elementary excitations from the solid to superfluid helium. 

Note, the resulting Hamiltonian contains the terms associated with creation or annihilation of single excitation as 
well as pair of excitations in helium, while it is only linear in creation and annihilation operators of phonons in the 
solid. This asymmetry is of course the shortcoming of our method. It can be avoided if we take into account the third 
order anharmonisms in the solid. However I am not concerned with this problem here because the main goal of the 
present paper is to describe the transmission of excitations from heated solid to the liquid helium at the temperatures 
closed to zero, which is consistent with real experimental situations. Under such conditions the contribution of single 
solid phonon annihilation is expected to dominate. 

III. TRANSFER HAMILTONIANS FOR PARTICULAR PROCESSES 

In this Section I write down the Hamiltonians for particular transfer processes: solid phonon O- single liquid phonon 
(single phonon emission), solid phonon two liquid phonons (two phonons emission), solid phonon single liquid 
roton (single roton emission). 

A. Single phonon emission 

To obtain the transfer Hamiltonian corresponding to the tunneling of a single phonon from the solid to the liquid 
helium and vice versa {i.e. to the processes q <-> k) we must retain in the expressions ( |l8| ) and ( |l9| ) only liner terms 
with respect to helium collective variable. In this case we get 



^ =Ti m + «vm n _^ (31) 



where I have used n = N/Ql, N is the total number of helium atoms. In the case of helium phonons = ck. Thus 
the relation (^2|) yields 

where I have neglected the term quadratic in wave vector. The same approximation can be done in the first term 
of eq. (|3"l|). In fact, without these approximations one can see that after substitution of eq. (^) into eq. © both 
small terms give even smaller one ~ (k 2 — k 2 ) 7 because k ~ k z in the process of interest. So, using the relations fl32] ) 
and ( pi] ) we come to 

Trf 1 / 11 ) =mc 2 ^n_ k e 4kr . (33) 

k 

Substituting eqs. (^5|) and d33| ) into eq. ( |l4| ) and integrating over the interface area S we obtain the transfer 
Hamiltonian for the single phonon channel in the following form 



G 



w£? h) =*Ecfc' PL 



1/2 f / , \ 1/2 



.k,q, V q,/ 



k,qt 



+ f^) V2 E *.(«■) f— ) 1/2 Mi. - 4U «4„.l. (34) 



k,q s 



This formula will be used in Section IV for deriving the Khalatnikov formula for heat flow via the single phonon 
channel. 



B. Two phonons emission 



The obtained formulae (|14|), (JlSf), and ( J19| ) allow us to study also the processes of emission of two helium phonons 
from single phonon in the solid. For this purpose let us retain in eqs. ( |l8| ) and ( |l9| ) only the terms quadratic in 
collective variables of liquid helium. This leads to the following expression for stress tensor 

= E [pLVteW, + g^n_ k n_ k j e« k + k > (35) 
which gives the transfer Hamiltonian for the two phonons processes under consideration in the form 

+ (^) 12 >. ^(,y,ll=^) E(k.k'.q,)S. C-!(i! 



£ B s ( 9s )(^) 1/2 5(k,k',q s )}. 

C:k',q 3 \ q s / J 



v ' k,k',q 3 

Here I have introduced the functions 

E(k,k',q a ) = Ji(q ,k|k / )<S q „ || +k||,k{ | - J^qal^k'^q^.kii+kj, (37) 

where 

1 k k ~ 



J 1 (q a ,k|k')= U + ^f) (^«k«L+^4«k'), (38) 



J 2 (q Q |k,k') = i + (^fitflt, + 6^,) . (39) 

Formulae ( |36| ) - ( |39| ) provide a possibility to study the energy transfer through the solid - liquid helium interface 
without conservation of number of quasiparticles. In this case there is no strong restriction on the angles for created 
phonons. Therefore such processes can-.lead to the isotropic background, which has been observed in experiments on 
phonon emission from cleaved crystalsBlii This problem will be considered in detail in another paper. 



C. Single roton emission 



To obtain the tunnel Hamiltonian, which corresponds to creation of a single roton in He II by a solid phonon, we 
can again use the relation ([H]) with account of eq. (22). In this case we however cannot neglect the terms quadratic 
in wave vectors. Note, the frequency w k in (|3l|) represents now the roton spectrum, whose analytic form depends on 
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the chosen concrete model. Here I do not distinguish the R^ - ) and R( + ' rotons. This problem will be considered in 
detail in Section V. 

Substituting eq. (E2) into eq. (J3l|) we get 



2^1.2 

ikr 



zz K ) L>\ k 2 Am J 



n_ k e lkr . (40) 



k 

Finally the corresponding transfer Hamiltonian takes the form 

1/2 



ki|,q, 



k,q t v ^ tJ 



2_, K ^Jl/2 (** b l - < b «>) *« •> ]• ( 41 ) 



where 



In Section V I will use the Hamiltonian (|4l]) for derivation of energy flux through the interface via the single roton 
channel. 

As the curtain fell of this section it should be noted that using the general formulae (|l8j) and we could study 
also various inelastic processes with rotons, e.g. solid phonon <^ two rotons, solid phonon <=> roton + phonon in 
helium, and scattering of roton on the surface with creation of a phonon in the solid. 

IV. SINGLE PHONON EMISSION. DERIVATION OF THE KHALATNIKOV FORMULA 

Any new theory claiming to describe the emission of elementary excitations from the solid to the superfluid helium 
should first be tested on the basic result obtained by Khalatnikov for a heat flow via the single phonon channeltx 
Really, this latter result was derived within appropriate framework, and is in excellent agreement with the experimental 
observations!™. 

To calculate the energy flux through the solid-liquid interface via single phonon channel we must first calculate the 
probability density function defined by the usual relation 



^(i P h) = 2hl 
m 2 



q<->k 



2 



<y(w q -Wk), (43) 



where — (q|7^^^ >l1 "' jk) is the corresponding matrix element of the transfer Hamiltonian given by eq. (|3J 

Substituting (|34|) into (f43|) we easily get 



(2tt) 3 {^ (lph) ,5(k|| -q i|| )+ ^ lpll) ^(k|| " q t || )+ 4 lph) <5(k|| -q s )}% q - Wk ), (44) 



where the coefficients A a are defined by the relations 

A (lph) = — -f^w. 4 lph) = Zr^rm, 4 lph) = ^^^fe 

To obtain formula (Q) I have done the transition from discrete to continuous wave vectors: 

<Sk||, qoi ,| — > ^- s ( k \\ - <M )■ 
The probability of creation of phonon with group velocity 
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"T5 («> 

near the surface element || e z (e z is the unit vector along zth direction) is given by 

^ (lph) ^ k n(Tg) (c ' dS) « M^ (lph) ?ia; k n(Ts)dS, (46) 
c 

where n(Tg) — [exp(?kj qcv / k^Tg) — 1] 1 is the equilibrium distribution function of phonons in the solid, Tg is the 
temperature of the solid. Thus, the heat flux from solid to liquid helium per unit area can be written as follows 

Q L d 3 k Q s d 3 qi 



f ^2 d 

f J A t fiuj ]s _n(Ts)5{^\\ -q t ||)(5(w q4 ~ w k)-^; 

y Aftw k n(T s )(5(k|| -q s )(5(w qs - o; k ) ^ -^-^ S 



) 3 (2tt) 3 

fJ L d 3 fc n s d 3 q t 
(2tt) 3 (2tt) 3 

(47) 



Delta functions of the type J(k|| — q Q ) take away integrations over k||. Delta functions S(ojq a — w k ) can be excluded 
by the integration over k z (note, k z ss k). After these straightforward calculations the integrals over the energies and 
over angles can be separated, so that the final result takes the form 

oo 

(i ph) A-kTi p L cF(rf) f u 3 duj n 2 p L c 4 

Qs^l = 7^3 §- / TT^ = onla 3 ( k b T s) F W> ( 48 ) 

(2tt) 3 ps s 3 J exp n^Jj _ ! 30 h 3 Ps s 3 

where F(rf) is the function of rj = si/st written down in Ref. 4. The formula ( fi"8|) exactly coincides with the well 
known Khalatnikov resultQ derived initially in the framework of acoustic mismatch theory. Note, the total heat flux 
between the solid and liquid helium is given by the difference 



j(iph) =Q (ipM_g(^) (49 ) 



where Q^E^s can ^ e f° un d in the similar way. Finally it is given by the same formula as (48) up to the replacement 
Ts — > Tjj , where Tl is the temperature of liquid helium. 



V. ROTONS EMISSION 



The novel techniques developed during last decade provide a possibility to produce and unambiguously detect by 
quantum evaporation processes both the R( — J rotons (k < fco) and the R^ + ) rotons (k > ho). However, it was realized 
that simple flat metal films produce easily measurable evaporation signals due to R( + ) rotonsc2l, but no signal from 
R'~) rotopSr-has ever been identified from such sources. To produce R^ - ) rotons two types of sources have been 
develop edcifE3 , which are based on trapping of the emitted R^ + ) rotons to enable R(' xotons to be created by 
interactions in nonequilibrium gas of excitations. Furthermore, Tucker and Wyatt reportedcj about creation of fast 
enough pulsed source of R^ - ) rotons suitable for time of flight measurements. 

Let us at first demonstrate the impossibility of direct emission of R' _ ) rotons from a solid, and, secondly, to describe 
the process of R( + ) rotons emission. 



A. R l ' rotons 



It is easy to see that R'~) rotons and R( + ) rotons arc not principally distinguished neither in the stress tensor 
( fl(]| ) nor in the corresponding transfer Hamiltonian (fill). Really the probability of R( _ ) roton emission seems to be 
nonzero if the energy and parallel component of momentum are conserved. This is just the case in the emission of 
real particles, say quantum evaporation^!. But quasiparticles in media are the collective excitations, which originated 
from local spatial perturbations. In this case the conservation laws at the surface should be obtained directly from the 
boundary conditions. In fact, the boundary conditions (0) and (R) just led us to the interaction Hamiltonian (0) we 



9 



have used. But the reversed conditions (which are clearly absurd) will necessarily give the same expression (||)! This 
means that for accurate use of this interaction Hamiltonian we must always keep in mind at least one of the boundary 
conditions, say (0) (the-second condition (^) is automatically satisfied as long as we use the transfer Hamiltonian (§])). 
In original SBT thcorjO there were no boundary conditions at all. This is the reason why the problem of R/ - ) rotons 
creation even did not appear in Rcf. 13. 

Let me now analyze the first boundary condition (^J). In view of the expression (|24| ) the normal component of 
displacement Sr z (x,y,0) at the surface can be written as followscil 

fr -<*''°> - (d^IX" ¥ <=""". m 

while retaining for simplicity one polarization, say longitudinal, eq. ( p5| ) at the interface can be rewritten in the 
analogous form: 

"•<^ > - (jSta)'"?^? ( S ' " 6L ') e ' ,,R "' (51) 



The right hand sides of eqs. (50) and ( pl| ) can be equal only if all their matrix elements in arbitrary basis are equal. 
Just this requirement leads to the conservation of the component of momentum parallel to the interface. But this can 
only occur if the coefficients at exponents in eqs. ( p0| ) and ( |5l| ) have the same sign! 

The detectable quasiparticle propagating into He II must have the positive zth component of group velocity. As is 
known momentum of R/ - -* rotons is oppositely directed to their group velocity. This results in impossibility to satisfy 
the boundary condition (||) in view of eqs. © and (g). 

B. R (+) rotons emission 

Let us now concentrate our attention in the processes of R( + ) rotons emission. In analogy with the previous section 
we can calculate the probability density function for single roton creation as 



W {1K) 



2tt 



2 



S(ljJq - W k ) 



= (2tt) 3 {A (1R) <Kk|| - eu||) + 4 XR) tQt\\ - qt||) + 4 1R) <5(k|| - q s )} «(w q - Wk), (52) 

where k denotes now the roton wave vector, co>k is its energy. The coefficients A\ are determined by the relations 

A (m)_ PL fc 2 B?(e)F*(k) , m) = PL fc 2 B?(6)F*(k) ^ 
1 ps^L^s 4m 2 w k w qi ' * psSIl^s 4m 2 w k w qt 

A{m = PL fc 2 5 2 fa)F 2 (k) _ 
The probability of creation of roton with energy ?Lwk and group velocity 

near the surface element || e z is equal to 

iy( 1R )^ k n(T g ) (Vr ' rfS) = H/ (1R) ^ k n(Ts)cos^E, (53) 

v r 

where cos 9 = v r • e 2 . Therefore, the energy flux per unit area associated with single roton channel takes the form 
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Q { ™{ = (2^) 3 | J 4 1R) ^ k cosen(T s )5(k|| - qi ||)5K, - Wk )^ 



fl L d 3 k Q s d 3 qi 

f (2tt) 3 

ft L d 3 fc fl s d 3 q t 



ytj 1R) 7kj k cos6>n(Ts)(5(k|| - q t ||)% qt -Wk)^ 3 ^ , a 
+ 1 4 1R Wos0n(T s )5(k|| - q.)«K "^T^f * "''^ ' 



r) 3 (2tt) 2 



(54) 



In principle the partial integrations in eq. (54) can give the spectral and angular distributions of emitted rotons. Here 
I calculate the total energy yield from single roton channel. The account of all phonon polarizations in the solid can 
only give the multiplier similar to F (rj) appeared in eq. (ff8|). This multiplier is of the order of unity for majority 
of solids. The account of all polarizations can have an essential influence only on the angular distribution. So, to 
calculate the total energy flux let us restrict ourselves to one polarization, say longitudinal one. Further I drop all 
the polarization subscripts. Thus, the energy flux of interest looks 

,(1R) TiPL f fc 2 ^ q ^k 



hpL f k 2 u>q0j. 2 cos 9 J mui 2 . Ti 2 (k 2 — k 2 



4(2Tr) 3 m 2 p s J exp ^^^-^ - 1 \ k2 4m 

X5(q|| - k\\)S(u}q - uj k )n(T s )d 3 kd 3 q. (55) 

After the straightforward, but quite cumbersome calculations eq. ( |55| ) can be reduced to the single integral. Let me 
write down here the spectral distribution for rotons emission: 



dk 120(2^)W S 7 P s exp (|^) - 1 t K 2 k* Vfc 4 < 



To obtain the total energy flux via the single roton channel Qs^l > ^ ne relation (^) should be integrated over k in the 



range {fco,/c max }, where fc max is the maximum roton wave vector. Note again that the formula (56) gives the energy 
flux in the case Tl = 0. In general, Ql^g can be calculated in the same way. 

To make a numerical analysis of the obtained results we must choose the particular analytic expression for roton 
spectrum in He II. Precise numerical calculations require of course the analytic expression most closed to a reality. 
But to draw the general picture I will use here the classical Landau expression^ 

h 2 (k — fc ) 2 /-»\ 
huj k = A+ v - u; 57 
2/i 

where A « 8.6 K, fc w 1.9 • 10 8 cm" 1 , (jt w 0.13m Hc 4 w 0.868 • 10" 24 g. 

Figure |^ presents the spectral distribution of roton emission calculated from (|56|) for gold. Figure clearly shows 
that the maximum at k > fco appears when temperature of the solid exceed 2 K. Integrating eq. (56) over R(+) 
rotons momenta we get the total energy flux. The ratio of this flux associated with the single roton processes and 
that corresponding to the single phonon processes given by eq. ( p8| ) is shown in Figure || as a function of temperature 
of the solid (again gold). Previously SBTE2I had also estimated the contribution of single roton emission processes 
into the total energy transmission through the solid - liquid helium interface. They, however, considered the case 
of small differences between T$ and Tl. Corresponding formulae can be easily obtained from eqs. ( [56] ) and ( pq ) by 
differentiation with respect to a temperature and subsequent integration over roton and phonon momenta, respectively. 
The corresponding ratio is shown in Figure [3] together with SBT result. The estimation by SBTE3 gives the ratio 
approximately three times larger for Tg = 2 K. This difference is caused by the approximations done in Ref. 13. 
Particularly the authors had taken k z w k, which is far from reality for rotons. Additionally they included Pj - ) 
rotons in their calculations, which is incorrect as I show in the previous subsection. But one can see from Fig. [| that 
their rough estimation nevertheless show right qualitative behaviour of the energy flux due to roton emission. 

VI. CONCLUSIONS 

I studied the emission of elementary excitations into liquid helium II from— a heated solid. For this process I 
developed the new microscopic quantum theory generalizing the previous theoryij. To describe the coupling between 



11 



the solid and He II the correct expression for the interface Hamiltonian (g) has been derived. This expression follows 
from the general principles of quantum mechanics with account of boundary conditions at solid-liquid, interface @j 
and (fj]). Following the ideas of Ref. 13 and using the microscopic description of superfluid heliumt3 the transfer 
Hamiltonian (^) has been expressed in terms of creation and annihilation operators of elementary excitations in the 
solid and helium II. The transfer Hamiltonians are explicitly written down for i) single phonon emission (|34|), ii) two 
phonons emission (^6|), and iii) single roton emission ([ill). The heat flux trough the interface associated with the 
single phonon channel is calculated by the use of the derived microscopic expression (|54[). The obtained formula ( |48| ) 
exactly coincides with the corresponding result obtained by Khalatnikovu using the acoustic-mismatch approach. The 
heat flux associated with the single R( + ) roton channel (see formulas ( p5| ) and ([56])) has been calculated as well. It is 
shown that the boundary conditions @ and (0) forbid the emission of R(~) rotons by a solid. This latter result is in 
accordance with experimental observations. 
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FIG. 1. Geometry of the problem 
FIG. 2. Spectral distribution of roton emission into He II from golden heater calculated by eq. fl56|) 



FIG. 3. Ratio of energy flux associated with the single roton channel and that corresponding to the single phonon (Khalat- 
nikov) channel in case Tl = 

FIG. 4. Ratio of energy flux associated with the single rotcp-ichannel and that corresponding to the single phonon (Khalat- 
nikov) channel in case Tl ~ Ts- Corresponding result of SBTtij is presented by the broken curve 
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